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INTRODUCTION 

A ballistic vehicle re-entering the earth's atmosphere at high speeds 
possesses a large amount of kinetic energy which is converted into aerodynamic 
heat. The terrific heat which surrounds the missile body requires that the body 
be coated with some type of material that can withstand this intense heat. Theo- 
retical and experimental investigations have proved that certain type glasses 
make a very good protective covering for a re-entry vehicle. 

For the approximately spherical surface in the vicinity of the forward stag- 
nation point of an ICBM, this report presents calculated results pertaining to 
the transient heating and ablation processes of 26 homogeneous, opaque, and non- 
decomposing glass shields, each shield having a unique set of thermal properties. 
The effect of the thermal properties on the heating and ablation results was 
investigated by varying each thermal property over a practical range of values. 

A similar study of the effect of the thermal properties on the heating and abla- 
tion results has been made for an IRBM re-entry [1, 4]. 

The results presented in this report for the ICBM re-entry were obtained by 
employing the numerical calculation method derived in references 3 or 5 and were 
evaluated on an IBM 704 computer. 


THE TRAJECTORY AND BODY CHARACTERISTICS 

The ICBM is assumed to begin its re-entry into the earth's atmosphere at an 
altitude of 150 km at Mach 15.68 or a flight speed of 7225.9 m/sec. Figure 2 
shows the flight altitude and flight speed as a function of time for the ICBM 
re-entry. The time t = 0 sec corresponds to the altitude of 150 km. The re-entry 
angle (0 e) is 114°, measured from the vertical, and the angle of attack (a) that 
remains constant throughout the re-entry is 0°. The horizontal range of the 
missile from zero time (H = 150 km) to impact time (H = 0 km) is 330.3 km. The 
ballistic factor (c fi ) of the ICBM is 1000 lb/ft 2 or 0.002 m 3 /kgsec 2 according to 
different definitions . 

The spherical portion of the glass shield at the stagnation point is assumed 
to have a diameter of 0.635 m. 


MATERIAL PROPERTIES AND FUNCTIONS EMPLOYED 


Since there is a lack of experimental data at high temperatures for the 
thermal conductivity (k) , the specific heat (c p ), and the emissivity constant 
(€), these material properties were treated as non-temperature-dependent con- 
stants for each calculated glass shield. The effect of the variation of these 
material property constants is, however, investigated by the calculation of 26 
supposedly opaque glass shields (see Table 1) with each glass shield having a 
different set of material properties. The specific weight ( 7 ) is assumed to be 
2105 kg/m 3 for all the calculated glass shields in this report. 

Since experimental data is available for the viscosity and the vapor pressure 
as a function of the temperature, these two material properties of glass were 
treated as temperature-dependent functions for each calculated glass shield. The 
viscosity functions for Pyrex (Corning 7900) and fused silica (quartz) are curve 
fits of curves presented in Figure 9 of reference 11. The vapor pressure function 
and heat of vaporization (h v ) for Pyrex were taken from reference 7. The vapor 
pressure function and heat of vaporization for fused silica were taken from 
reference 9. When the constants p* and A in the viscosity and vapor pressure 
functions are equal to one, the equations given below for p and p v are the exact 
curve fits of these two properties for Pyrex and fused silica glasses. 

The vapor pressure function, viscosity function, and heat of vaporization 
for Pyrex and fused silica glass are as follows: 

(1) Pyrex glass 

(a) the vapor pressure function 


Pv, P = 10>332 A p exp 
(b) the viscosity function 

Pp = 0.016684 p'p exp 


-46,400 , 

T(°K) + 14 ' 5 


4230 

T(°K) 


2.0635 


(c) the heat of vaporization 
h V)P = 2470 (kcal/kg) 

(2) Fused silica (quartz) 

(a) the vapor pressure function 

-57,800 


(kg/m 2 ) 


(kg sec/m 2 ) 


P V ,F.S. = 10 > 332 ex P 


T( K) 


+ 18.48 


(kg/m 2 ) 


( 1 ) 


( 2 ) 


(3) 


2 



(b) the viscosity function 


Up s = 4.9505 exp 


4836. 97\ 2 * 6366 

. T(°K) J 


(kg sec/m 2 ) 


(4) 


(c) the heat of vaporization 
h v,F.S. = 3050 ( kcal / k §) 

The constants Pp and Ap are one for the real material and were varied in some of 
the cases studied here. 

The material properties that were varied and the limits within which each 
material property is varied in this report are as follows: 

(1) the thermal conductivity, k (kcal/m °K sec) 

0.61 x 10 -4 < k < 48.84 x 10~ 4 

(2) the specific heat, c p (kcal/kg °K) 

0.10 < c p < 0.58 

(3) the emissivity constant, e (-) 

0.40 < e < 1.0 

(4) the vi scosity level, (ip (-) 

0.10 < |Jp < 00 

(5) the vapor pressure level A p 
0 < A p < 1000 

Since the cases as the individual material properties k, c p , and p* approach zero 
and infinity were discussed in reference 1, they will not be repeated in this 
report. 

A so-called '’standard set 11 of material properties will be referred to in this 
report, and they are as follows: 

k = 0.61 x 10" 4 (kcal/m°K sec) 

Cp = 0.29 (kcal/kg °K) 
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e = 0.8 

(-) 

- i.o 

(-) 

II 

o 

(-) 


The "standard set" of material properties was employed in the theoretical study 
of the ablation of a glass shield at the stagnation point of a re-entering IRBM 
[l] and a re-entering satellite [6], 


THE EFFECT OF THE MATERIAL PROPERTIES ON THE INTERFACE TEMPERATURE 

The interface temperature T^ follows the same trend with respect to time for 
each glass shield regardless of the material properties, and the maximum inter- 
face temperature occurs at approximately the same time for each shield since 
h e » h Therefore, the effect of the material properties on the interface 
temperature can be studied by the consideration of the maximum interface tempera- 
ture (T^ max ) for the ICBM re-entry. Figure 3 shows the effect of the material 
properties on the maximum interface temperature. Large values for the thermal 
conductivity k, specific heat Cp , emissivity constant e, and vapor pressure level 
A tend to keep the surface temperature low; the viscosity level |i* has very little 
effect on the surface temperature. 

Figure 4 is a pictorial view of a Pyrex glass shield (Case 24, Table 1) 
showing its complete temperature and ablation history for the ICBM re-entry. One 
can see the movement of the interface at the melting rate -v GO (ram/sec) toward 
the thermally insulated inner edge of the glass shield in this graph. 


THE HEAT FLUXES 
1. ilie Aerodynamic Heat Flux 

The compression at and near the forward surface and the friction between 
the fluid particles as they flow along the body converts the kinetic energy of 
motion into aerodynamic heat. The formulas that were used for the aerodynamic 
heat flux rate at the stagnation point of a non-ablating surface and the source 
of these formulas are as follows: 


(1) Molecular flow region [10] 


q = 1.054 x 10" 3 p V 
aero «■> oo 


oo 


V" - 1.9263 x 10 3 (T. - T ) 

1 oo 


kcal/m 2 sec (5) 
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(2) Continuum flow region for V > 2100 m/sec [8] 


‘aero 


h - h. 


q = 74,571.6 n/poo/r' (V / u) 3,15 (t ^ 77 ) kcal/m J 

9 00 c V h 73/ 

v e 

(3) For V < 2100 m/sec [12] 


sec 


( 6 ) 


oo 


q = OcA/TT) • • (T - T.) -s/T p /T.' (Nu /Re7) D , (.715)°’ 4 -/iC. • V /D* 

aero v 7 ^ ; i,air e 1 / e'e l v 1 Pr = 1 M co 


kcal/m 2 sec , 
(7) 


The aerodynamic heat flux for an ablating surface, e.g. , glass vaporization, is 
given by 


q = q \j r 

aero aero 


where ^ is given in Reference 7 as 


t - 


1 - c 


1 - (1 - Kj. M ni )c 


( 8 ) 


(9) 


c 


1 

1 + M (2s- - 1) 

Pv 


( 10 ) 


For the glass shields considered herein, M = 0,72, K x = 0.68, and n 1 = 0.26, 
The relation for the surface ablation taken from reference 7 is 


(7v). = 


- q c 

^aero 


1 (h e - h.) [1 - (1 - K x M ni )c] 


( 11 ) 


The time t = c (sec) designates the time after ablation has ended when the 
net heat flux 


q . « q -q , + h ( 7 v) . ( 12 ) 

l aero rad v ' 'i 

across the surface changes its sign from positive to negative. Beginning then, 
the radiative heat transfer q rac | out of the shield exceeds the aerodynamic heating 
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rate T aero * The time t" is between the time of peak heating and impact time and 
is approximately the same for all the glass shields. All of the individual heat 
fluxes for each glass shield were integrated over time t to t = t'\ since 

t* 

Q t (t*) = J q. dt 

o 

gives the maximum amount of heat that the shield, in both the solid and the liquid 
stage, will absorb. The time integrated q f s to t* will be designated Q(t*) in 
this report. 


The aerodynamic heat flux Q aero (^ C ) that stands for the time integrated 
aerodynamic heat flux without mass flow across the interface is given for each 
one of the calculated glass shields in Table 1, The true aerodynamic heat flux, 
QaeroV*)’ is given by 


Q (t*) = Q (t*) 
^aero aero 


Q B1 + 


h (-/v) . dt, 
v ' 1 


(13) 


where Q B l(t') is the total amount of aerodynamic heat blocked up to time t' c due 
to the vaporization process of the glass. Table 1 shows that the total aero- 
dynamic heat, Qaero(t*)> changes very little for all of the calculated cases 
since h e » hi during the period of significant heating. Thus, the material 
properties of the glass shield have very little effect on the aerodynamic heat 
flux Qaero* The trajectory of the re-entering vehicle principally determines 
the aerodynamic heat flux Q aero * 


2. The Heat Flux Radiated Away From the Opaque Glass Shields 


The heat flux that is radiated away from the outer surface of the sup- 
posedly opaque glass shield given by 


q j = 1.378 x 10" 11 € T? (kcal/m 2 sec) 
rad l ' 


(14) 


is a function of the wall temperature and emissivity constant e. The values of 
Qrad(t^) listed in Table 1 show that the outer surface radiates away the major 
portion of the aerodynamic heat flux Q a ero( tX ) for most °f the calculated shields. 
The radiative efficiency of the glass shields is defined as the fraction 
Q r ad( t ) /Qaero (t*) of the gross heat input that is radiated away from the glass 
wall. Figure 5 shows the radiative efficiency of the heat protection shields for 
the varied material properties assumed for the glass. This graph shows that the 
radiative efficiency of a glass heat shield is optimized by utilizing a glass 
with a small thermal conductivity and specific heat and with a small vapor pressure 
Changes in the viscosity level and emissivity constant have a much smaller effect 
on the radiative efficiency than the other pertinent material problems. 
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3. The Heat Blocked, Qg^(t^), Due to the Vaporization Process of the Glass 


The total amount of heat blocked due to the vaporization is given by 


<B1 


= Q 


aero 


aero 


+ h 


(-/v) i dt, 


(15) 


where (Qaero “ Qaero) represents the amount of heat blocked by the diffusion of 
the vapor across the boundary layer and 

t" 

h v J (-rv) i dt 

o 


is the amount of heat absorbed by the evaporation process at the surface of the 
shield. Cases calculated for Ap = 1 , 3 , 100, and 1000 (Table 1), with all other 
material properties being the so-called "standard set" of material properties, 
resulted in 46 to 88 percent of the total aerodynamic heat Q aero ( t *) being blocked. 
The amount of the total ablation that was due to vaporization ranged from 33 to 
98 percent. Reference 5 points out that radiative cooling is more desirable than 
mass transfer cooling under certain conditions. 


4. The Heat Flux Qp(t'') Absorbed by the Glass Shield 


The total amount of heat that the molten and solid layers of the shield 
absorb is given by 


■ <Lro <t * ) ' Q rad (t *> ‘ V (t *> " 


■ Q (t*) - Q ,(t*) - h / (-yv). dt. 
aero rad v / i 


(16) 


This total heat absorbed Qj_(t*) per unit area by the shield ranges from 3146 kcal/m 2 
to 60^000 kcal/m 2 for all cases calculated. The amount of the total aerodynamic 
heat Qaero (t*) absorbed by the shield ranged from 7 to 91 percent. 

Figure 6, presenting Qj_(t*) as a function of the material properties, shows 
that the thermal conductivity and specific heat affect the magnitude of absorbed 
heat Qp(t*) by the shield much more than the other material properties. The 
amount of heat Q-^(t x ) should be small since then the shield thickness needed as 
a heat sink is kept to a minimum. Q^(t' v ) can be kept small by employing a glass 
with a small thermal conductivity, a small specific heat, a high viscosity and 
vapor pressure level, and a large emissivity at the surface. 
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5. The Amount Q c (t*) of the Absorbed Heat Carried Away by the Molten Flow 

The convective heat flux Q c (t ) is the amount of heat that the flow of 
molten glass carries away from the stagnation point vicinity. This heat Q c (t*) 
carried away by the liquid glass is part of the heat Qi(t*) that. was absorbed by 
the glass shield. Table 1 shows that from 0 to 79 percent of the total aero- 
dynamic heat Qaero^*) carried away by the molten flow. The amount of the 
heat absorbed Q^(t*) that is carried away by the molten flow ranged from 0 to 
96 percent. Table 1 and Figure 6 show that those shields which absorbed a large 
percentage of the total heat QaeroC* 1 *) h a< 3 most of this absorbed heat Q-j_(t*) 
carried away by the molten flow. 

The difference - Q c (t*) is that amount of heat the glass shield 

will have to store in heat sink fashion at the stagnation point, because after 
time t' no more convective cooling takes place since melting has stopped. After 
time t* the heat Qi(t*) - Q c (t*) is very important since this heat keeps pene- 
trating into the shield until impact time tf. The fraction R(t*) given by 

(Qi(t*) - Qc(i*)) /Qaero^*) is the fract i° n of the gross heat Qaero^*) that is 
stored in heat sink fashion in the shield. This fraction R(t*) varies from 0.027 
to 0.122; this shows that the glass shields have to store very little of the gross 
heat . 


THE ABLATION OF THE GLASS SHIELDS 

The ablation rate, -v 00 (mm/sec), of the glass shield is due to (1) the 
flow of molten glass and (2) the evaporation of the glass at the surface. The 
total ablation s, i.e., the layer lost due to ablation, given by 



for each glass shield, is given in Table 1. For the shields investigated, the 
thickness s of the heat shield lost due to ablation was from 0 ,to 35 mm, amounting 
to less than 1.4 inches. 

Figure 7 shows the total ablation s as a function of the material properties 
of the glass shield. Figures 7, 8, and 9 show that the thermal conductivity k, 
and the specific heat Cp are the material properties that principally determine 
the total ablation s when the ratio k/yCp is small, as is desirable. Practical 
variations of the viscosity level pi*, surface emissivity e, and the vapor pressure 
level A caused rather small variations in the total ablation for the glass shields 
which had a small thermal diffusivity k/ycp. 

The dependency of s on the viscosity level for a non-vaporizing glass 
shield is shown in Figure 8 by comparison of s for cases 8, 19, and 20. All three 
cases have the properties k = 9.7 x 10“ 4 kcal/m °K sec, Cp = 0.29 kcal/kg °K, 
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k/ycp - 15.9 x 10“ 7 m 2 / sec , and A p = 0. The viscosity level |i* p was varied within 
the limits 0.1012 < |i*p < 101.2 which resulted in the total ablation s ranging 
between 9.3 < s < 35.7 ram. Evidently, s increases with a decreasing viscosity 
level. 


The effect of varying k and c p for a large constant ratio k/ycp is shown 

by comparing Cases 21 and 22 (Fig. 8) where k/^Cp = 40 x 10” 7 m 2 / sec for both 

cases; Case 21 has values of k and c p which are smaller than those of Case 22 

by a factor of two. Case 21 with the smaller k and c p resulted in an ablated 

layer 59% larger than that of Case 22. 

Cases 2, 10, 11, 25, and 26 in Table 1 and Figure 7 show the effect on the 
total ablation s of varying the vapor pressure level Ap between the limits 
0 < A p < 1000. The total ablation s increases only slightly with increasing 
vapor pressure level. Table 1 shows that the fraction 



of the total ablation due to vaporization is as large as 98% for Ap = 1000, and 
that this fraction increases together with the vapor pressure level A and vis- 
cosity level |i*. 


THE THERMAL PENETRATION INTO THE GLASS SHIELD 

The calculated "semi-infinite solution" [3] is based on the condition that 

lim T(t, z) * T 0 
z — > 

where T 0 was assumed to be 300°K in this report. For all the cases considered, 
a thermally insulated inner wall of the shield was assumed to have been placed 
at the maximum distance (denoted z = b) from the final surface where the semi- 
infinite solution yielded T = 380°K. For the given re-entry, the time of maxi- 
mum thermal penetration was impact time tf. The true temperature at the point 
b and time tf is 

T( t f , b) = T 0 + 2(380 - T 0 ) = 460% (18) 

which follows from adding the temperature profile given by the semi-infinite 
solution to the reflection of this profile about the point b. This satisfies 
the conditions that &T(tf, b)/5z = 0 and that the temperature of the thermally 
insulated inner wall is < 460°K. The numerical results show that this assumption 
causes a negligible temperature increase at the outer surface. 
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The maximum thermal penetration b, or the thickness of glass needed as a 
heat sink, is given in Table 1 for all of the calculated cases and is plotted as 
a function of the material properties in Figure 7. The thickness b varies from 
2.6 to 21.15 mm (Table 1). Figure 7 shows that the thermal penetration b is 
predominantly governed by the thermal conductivity k and specific heat c p . The 
viscosity level J r has some effect on b, but not as much as k or c p , whereas the 
emissivity e .of the surface and the vapor pressure level A had an almost negli- 
gible effect on this thickness b (Fig. 7 and 8) since large percentages of the 
total radiation emitted from the surface and the total transport of material by 
melt flow and evaporation are concentrated in the short -lasting, high-temperature 
range of the surface. In this range, small changes in the surface temperature 
can compensate for large changes in the emissivity constant and in the assumed 
level factors A and |j* of the vapor pressure and viscosity, respectively, because 
of the highly nonlinear relations p v = p v (T) , \i = fl (T) , and q rad = q rad (T) . 

Figure 8 shows the thermal penetration b as a function of the thermal 
diffusivity k/y c p for all the calculated glass shields; it is seen that b 
increases together with k/7c p and fi" 

THE NECESSARY WEIGHT OF THE GLASS SHIELDS 

The necessary weight per unit area of the glass shield, defined as y ( s + b) , 
is the most important parameter of the ablation type heat shield to be considered; 
the sum (s 4- b) is defined as the necessary thickness. The weight of the protec- 
tive glass shield should be kept as light as possible, but the shield should 
still successfully protect the re-entry body. The specific weight y of the glass 
shields was assumed to be constant, i.e., y = 2105 kg/m 3 , in this report; 
therefore, the two constituents that determine the necessary weight y (s 4- b) are 
the total ablation s and the heat sink thickness b. 

For an IRBM re-entry [4], the thermal diffusivity k/yc p was the dominant 
material parameter that determined the necessary thickness (s 4 b) of a glass 
shield in the vicinity of the stagnation point. Melting and significant aero- 
dynamic heating stopped about 40 seconds before impact time for the IRBM re-entry. 
Since s was considerably smaller than b, the dependency of (s 4- b) on the material 
properties other than k/ 7 c p was small. 

For the ICBM re-entry at hand, melting and significant aerodynamic heating 
stopped only 15 seconds before impact time; also, the contribution of the total 
ablation s to the necessary thickness (s 4- b) is quite significant except for the 
cases where the viscosity level was assumed to be so high that little melting 
occurred. Therefore, (s 4- b) is affected by other parameters in addition to 
k/Xcp. Figure 9 presents the necessary weight as a function of the thermal dif- 
fusivity k/7c p for all the glass shields calculated. From the results it is 
concluded that the following trend in the material properties would provide the 
lightest glass heat protection shield for the ICBM re-entry: 
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(1) a small thermal diffusivity k/ycp , 

(2) a high viscosity level, i.e., a large resistance to melting, 

(3) a large emissivity e of the supposedly opaque glass. 

The results show that the vapor pressure level A has a negligible effect on the 
necessary weight y (s + b) . 


CONCLUSIONS 

The heat protection at the stagnation point of an ICBM re-entering the 
atmosphere of the earth is investigated in this report by using homogeneous, 
opaque glass shields. The effect of the material properties of the glass on the 
performance parameters is analyzed from the results of 26 calculated glass 
shields, each with a different set of assumed material properties. 

The most important performance parameter of the heat protection device is 
the necessary weight 7 (s + b) per unit area that is needed to keep the thermally 
insulated inner edge of the shield below an arbitrarily chosen low temperature 
of 460°K. The following trends of the glass thermal properties would tend to 
minimize the weight of the glass heat protection shield for the ICBM re-entry: 

(1) a small thermal diffusivity k/^Cp, 

(2) a high viscosity level, i.e., a large resistance to melting, 

(3) a large emissivity constant € of the supposedly opaque glass; 

this is not as important as the thermal diffusivity or the viscosity 
of the glass . 

It was found that the vaporization of the glass had very little effect on 
the necessary weight 7 (s + b) . 
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TABU 1. THE CALCULATED GLASS SHIELDS FOR THE ICBM 
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FIG. 1. THE REFERENCE SYSTEM 
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FIG. 4. TEMPERATURE AND ABLATION HISTORY OF A GLASS SHIELD 
AS A FUNCTION OF THE MATERIAL PROPERTIES 




FIG. 5. THE RADIATIVE EFFICIENCY OFTHE GLASS SHIELDS 
AS A FUNCTION OF THE MATERIAL PROPERTIES 
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FIG. 6. THE HEAT ABSORBED Qj ft* ) BY THE GLASS SHIELDS AND THE 
AMOUNT Q c (t*) OF THE ABSORBED HEAT THAT IS CARRIED 
AWAY BY THE FLOW OF THE LIQUID GLASS AS A FUNCTION 
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FIG. 8. THE TOTAL ABLATION s AND THERMAL PENETRATION b AS A FUNCTION OF 
THE THERMAL DIFFUSIVITY k/yc FOR ALL THE CALCULATED GLASS SHIELDS 







m 



023 


C 


OH 


s 


00 

Q 


o 

(SI 


o 


o 

~o 


o 

sO 


o 


23 

bi-e 

(0 

OQB0 

« 

o 

00 


10 

£L 


o o 

| CO 4 ) 
OB 


X 
(X 
u 

is 5 


o 

^ oo 

O °° 
— oo 

£3 

2 O 

£ Q 
uu 

<£ 
oo Zj 
<: zd 
_ o 

3 -i 
+ < 
^ O 

>-UJ 

S2=l 

uj C 

3 = a; 

se 
< >■ 
00 
oo 


o 


> 

00 


h- Q 

o 

o 


21 




REFERENCES 


1. Adams, E. W., "Theoretical Investigations of the Ablation of a Glass-type 
Heat Protection Shield of Varied Material Properties at the Stagnation 
Point of a Re-entering IRBM," NASA TN D-564, January 1961. 

2. Adams, E. W., "Radiation Versus Mass Transfer Effects for Ablation Re-entry 
Shields of a Nonlifting Satellite," Journal of the Aero/Space Sciences , 
Readers Forum, Vol . 27, No. 11, November 1960. 

3. Adams, E. W., "Calculation Method for the Transient Melting Process of a 
Glassy Heat Protection Shield in the Vicinity of the Stagnation Point of 
a Re-entering Vehicle," ABMA Report DA-TR-8-60, March 1960. 

4. Adams, E. W., "The Thickness of a Melting Ablation Type Heat Shield, 

Journal of the Aero/Space Sciences , Readers Forum, Vol. 27, No. 7, July 1960 

5. Adams, E. W., "Analysis of Quartz and Teflon Shields for a Particular 
Re-entry Mission," Proceedings of the 1961 Heat Transfer and Fluid 
Mechanics Institute , Stanford University Press, 1961, pp . 222-236. 

6. Adams, E. W., "Melting-Type Heat Protection of a Ballistic Satellite 
Re-entering the Atmosphere of the Earth," Fusees et Recherche Astronautique 
et Spatiale , Paris, pp . 7-20. 

7. Bethe, H. A., and Adams, M. C., "A Theory for the Ablation of Glassy 
Materials," Journal of the Aero/Space Sciences , Vol. 26, June 1959, 
pp. 321-328. 

8. Detra, R. W., Kemp, N. H., and Riddell, F. R., "Addendum to Heat Transfer 
to Satellite Vehicles Re-entering the Atmosphere," Jet Propulsion , 

pp .’ H256-1257 , December 1957. 

9. Georgiev, S., "Unsteady Ablation," Semi-Annual Report on the Re-entry 
Problem, AVCO Corporation, December 1958. 

10. Oppenheim, A. K., "Generalized Theory of Convective Heat Transfer in a 
Free Molecular Flow," Journal of the Aero/Space Sciences , Vol. 21, 

July 1954, pp. 459-474. 

11. "Properties of Selected Commercial Glasses," Corning Glass Works, Corning, 

New York, 1959 . 

12. Reshotko, E., and Cohen, C. B., "Heat Transfer at the Forward Stagnation 
Point of Blunt Bodies," NACA TN 3513, July 1955. 


22 


NASA- Langley. 1963 Ml03 



